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ABSTRACT 

We present the simplest model that permits a largely analytical exploration of the 
m = 1 counter-rotating instability in a "hot" nearly Keplerian disc of collisionless 
self-gravitating matter. The model consists of a two-component softened gravity disc, 
whose linear modes are analysed using WKB. The modes are slow in the sense that 
their (complex) frequency is smaller than the Keplerian orbital frequency by a factor 
which is of order the ratio of the disc mass to the mass of the central object. Very 
simple analytical expressions are derived for the precession frequencies and growth 
rates of local modes; it is shown that a nearly Keplerian disc must be unrealistically 
hot to avoid an overstability. Global modes are constructed for the case of zero net 
rotation. 
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1 INTRODUCTION 

Galactic nuclei are thought to harbour supermassive 
black holes and dense clusters of stars, whose struc- 
tural and kinematical properties appear to be corre- 
lated with global galaxy properties jG ebhardt et al.l fl996l : 
iFerrarese fc Merrittl l200d : iGebhardt et al.l l2000h ~The im- 
print of galaxy formation is expected to be recorded in 
the nature of stellar orbits. A remarkable case is that of 
our nearest large neighbouring galaxy M31, whose cen- 
tre has a double-peaked distribution of stars lILTght et al.l 
1974 lLauer et al.lll993l. Il998l : iKormendv fc Benderlll999T T 



19881; iMerritt fc Stiavehi 


1990; iPalmer & Pa 


paloizoulll990l; 


Sellwood fc Merritdll994; 


Lovelace et ah 1997 


). The dynam- 



ics of galactic nuclei involve nea rly Keplerian systems of 



stars or other collisionless matter dRauch fc Tremaine 



Sridhar fc Toumalll999l; ISridhar et al.lll999l ; lTremaine 



1990 



2001 



Touma et al.l 120091 ) ■ IToumal i^Ooj Tconsidered a softened 



Tremainel (| 19951) proposed that the off-centered peak marks 
the region, in a disc of stars, where lie the apoapses 
of many eccentric orbits. Self-gravitati ng models of such 
an eccentric disc have been proposed |Bacon et all l200ll; 
ISalow fc Statlerll200ll ; ISambhus fc Sridharll2002l ) . Of partic- 
ular inte rest to the present investi gation is the model pro- 
posed in ISambhus fc Sridhar! (|2002l ). because it included a 
few percent of stars on retrograde (i.e. counter-rotating) 
orbits. Here, it was proposed that the lopsidedness of the 
nuclear disc of M31 could have been excited by the presence 
of the retrograde stars, which were accreted to the centre 
of the galaxy in the form of a globular cluster that spiraled 
in due to dyna mical fri c tion. This proposal was motivated 
by the work of IToumal l|2002l) , which suggested that even 
a small fraction of mass in retrograde orbits could excite a 
linear lopsided instability. 

Counter-rotating streams of matter in a self- 
gravitating d isc are known to be unstab l e to lop- 
sided modes l|Zang fc Hohll 1 19781 ; lArakil 1 19871 ; ISawamural 



gravity version of Laplace-Lagrange theory of planetary mo- 
tions, and showed that a nearly Keplerian axisymmetric disc 
is linearly unstable to a m = 1 mode when even a small 
fraction of the disc m ass is counter- rotating. Softened grav- 
ity was introduced bv lMillerl { 1971) to simplify the analysis 
of the dynamics of stellar systems. In this form of interac- 
tion, the Newtonian l/d gravitational potential is replaced 
by 1/ \/d 2 + b' 2 , where b > is called the softening length. In 
the context of waves in discs, it is well-known that the soft- 
ening length mimics the epicyclic radius of stars on nearly 
circular orbits. Therefore, a disc composed of cold collision- 
less matter interacting via softened gravity provides a sur- 
rogate for a "hot" collisionless disc. 

The goal of this paper is to formulate and analyse the 
counter-rotating instability in the simplest model of a "hot" 
nearly Keplerian disc of collisionless self-gravitating mat- 
ter. To this end we make the following choices: (i) The discs 
are assumed to be made of matter whose self-interaction is 
through softened gravity; (ii) A Wentzel-Kramers-Brillouin 
(WKB) analysis is made of the linearised equations gov- 
erning the perturbations. The unperturbed two-component 
nearly Keplerian disc is introduced in § 2 and the apse 
precession rate is defined. The equations governing the lin- 
earised perturbations and relevant potential theory for soft- 
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ened gravity is given in § 3. The local (or WKB) approxima- 
tion and dispersion relation for local modes is derived in § 4. 
This is used to discuss stability, instability and overstability. 
§ 5 considers the construction of global modes for the case 
of a Kuzmin disc with equal masses in the counter-rotating 
components, and we conclude in § 6. 



2 THE UNPERTURBED TWO-COMPONENT 
DISC 

We consider a disc of mass Md orbiting a central mass M. We 
specialise to the nearly Keplerian case, Md <C M . Therefore 
the force on the disc material is mostly Newtonian, giving 
rise to a near equality between the frequencies of azimuthal 
and radial oscillations. Test particle orbits may be thought 
of as osculating Keplerian ellipses, whose apsides precess due 
to the self-gravity of the disc, at rates that are smaller than 
the orbital frequency by a factor ~ e — Md/M <C 1 . Other 
forces could also be responsible for the evolution of the disc 
over similar slow time scales. We consider a cold collisionless 
disc, composed of particles orbiting on two counter-rotating 
streams. The disc particles interact with each other through 
softened gravity. The central mass and the disc attract each 
other through ( unsoftened) New tonian gravity. Our notation 
closely follows iTremainel (|200ll ). 

We employ polar coordinates (r, (f>) in the disc plane and 
place the central mass at the origin. The unperturbed com- 
ponents are axisymmetric with surface densities Sj(r) ^ 
and E^ (r) ^ 0, and circular velocities = ±rfl(r)0, re- 
spectively. The angular frequency, fi(r) > , is determined 
by the total gravitational potential: 



$(r) 



GM 



+ $ d (r) + $ e (r) , 



(1) 



where <&d is the self gravity of the disc determined by the 
total surface density £<j(r) = E^J" + E^ , and $ e is a non- 
Keplerian potential due to an external source. <&d is 0(e) 
compared to GM/r; we assume that $ e is also 0(e) . Test 
particles on nearly circular prograde orbits have azimuthal 
and radial frequencies, Q > and k > ; particles on 
nearly circular retrograde orbits have frequencies, — fl < 
and — k < . The frequencies are given by, 



„2 1 \ GM Id,. 

fi 2 (r) = -j- + -— $ d + $e 
r J r dr 



GM 



d 2 



3 d 



(2) 



(3) 



The precession rate of the apsides of a nearly circular orbit 
of angular frequency ±fi(r) , is given by ±zu , where, 



zu(r) 



= Q, — K 
1 



d~ 



2 d 



2fi(r) \dr 2 r dr 



it, + :x (** + *«) + 0(0- 



(4) 



Euler and continuity equations, appropriate to a cold disc: 



On 
~dt 

at 



2- + (v$-V)v± + (i£-V)»J = -V$«, (5) 



(6) 



where & a (r,t) is the perturbing potential. We write 
the variables w*, v^, 3>a) in the form X a {r,cj),t) = 

ex P — wi)] ■ Substituting these in equations ((5]l 

and ©, straightforward manipulations give, 



u„ = — - 



(±m\l — uj)— ± 

dr r 



, k d m , , _ . 

± 2«^ + 7 (±mn -^ 



(7) 
(8) 



i(±mfi-a))S* 



r dr 



where 



Of 



(9) 



(10) 



Equations ((TJ — (|10[) give the linear responses of the surface 
densities and velocities of the two components, to a specified 
perturbing potential, "!> a . 

For a self-consistent response, the perturbing poten- 
tial, $™(r), depends only on the total surface density, 
(E™ + + E" l_ ); in fact this is the only coupling between 
the two counter-rotating components. The Poisson integral 



#™(r) = / dr' r' P m (r,r') [X? + (r') + E™"(r')] . 
Jo 

(11) 

The unperturbed disc potential, $d, is related to the un- 
perturbed surface density, Ed = (Ed + Ed) through equa- 
tion when m — 0. The kernel, 

P m (r,r') = -—BT(a,/3) + nG^(S ml + 5 m> _0 , 
r> r 

(12) 

includes direct and indirect contributions. Here r< = 
min(r, r'), r> = max(r, r'), a = r</r> , /3 — b/r > , and, 



B?{a,P) = 



d6 cos m# 



irjo (1 - 2a cos + a 2 + f3 2 ) 



2W 2 



(13) 



is a gener a lisatio n of the Laplace coefficients, introduced 
by iToumal (|2002l ): in the limit of no softening, B™(a,0) = 
b™ 2 ( a )> the Laplace coefficie nts familiar from celestial me- 
chanics (|Murrav fc Dermott|[l999l ). Equations (0 — 113| de- 
termine the self-consistent, linear modes of axisymmetric 
discs, whose counter-rotating components interact through 
softened gravity. 



3 LINEAR RESPONSE 



Let = Ua(r,t)r + v^{r,t)<j> and E„ (r,i) be infinitesi- 
mal perturbations toVj and E^. They satisfy the linearised 



4 THE LOCAL APPROXIMATION 

A perturbation, X?(r) = \X a n (r)\ exp [i f dr k(r)] , is re- 
ferred to as tightly-wound, if the radial wavenumber is large 
in the sense, [rfc(r)| 3> |m|. To leading order in \m/rk(r)\, 
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the WKB approximations to the linear responses of equa- 
tions 0— © are, 



u a (r) = 



(±mfi — uj) 

51 



k<Z>" 



= (r) 



e^M 



k^l 



(14) 



(15) 



(16) 



The responses are singular at radii, where = 0. For 
a self-consistent $™, the potential theory of th e previous 
Section simplifies in the WKB limit i|Millerlll97il ). 



exp (— \k\b ) 



(17) 



Substituting for E™ + and E™ from equation (|16[l . and 
eliminating $™, gives the WKB dispersion relation, 

D+Z?- = 2^G|fc|exp(-|fc|6) [D~S+ + 73+ E^ 



(18) 



which is a quartic equation in ui. All the coefficients of the 
various powers of ui being real, if a; is a solution, so is its 
complex conjugate, u>* . It should be noted that we have not 
made any assumptions about the Keplerian nature of the 
disc. Therefore the dispersion relation of equation (|18[1 is 
valid for a non Keplerian disc as well, with or without the 
central mass. In general, uj will not be small, compared to 
either Q or re, so the the dispersion relation is not restricted 
to stow perturbations. When the counter-rotating compo- 
nent is absent (i.e. E^" = 0), equation (|18p reduces to the 
WKB dispersion re l ation, familiar from Problem (6-5) of 
iBinnev fc Tremaind l| 19871 ). 

We recall the result for the stability of axisymmetric 
(m = 0) perturbations. When m = 0, we have D + — 
uj 2 , and equation (|18|l becomes, 

(19) 



27rGE d |fe| exp(-|fc|6) 



which is identical to the dispersion relation for a disc without 
counter-rotating co mponents. It is straightforwa rd to prove 
(as Problem (6-5) of IBinnev fc Tremainei i| 19871 ) invites the 
reader to) that the disc is stable to short-wavelength ax- 
isymmetric perturbations, if 

b>bo ^^£. (20 ) 



4.1 Slow m = 1 Perturbations 

When the azimuthal wavenumber is m = 1, the near equal- 
ity between f2 and re — see equations ([2} and ([3} — enables 
slow modes, for which \u)/fl\ — 0(e) <C 1. It is convenient 
to define, 

E+ = (l-r,)E d , and E^ = r?E d , (21) 

where ^ rj(r) ^ 1 is the local mass fraction in the counter- 
rotating component. We also introduce a local frequency, 



5(r,|*|) = *GE£)|*|exp(- 
whose maximum value at any r, 



c(r) = 



2Q,b 



(22) 



is attained for |fc| = 1/b. Here bo is the minimum softening 
length, defined in equation (|20p , that ensures local stability 
to axisymmetric perturbations. 

From equations (|4]) and (|10|) . we have, 

(24) 



DJ = 2Q (uj - vo) + Q(e 



-2f2(w + vj) + 0(e ). 



(25) 



When these are substituted in equation (|18p . a little rear- 
rangement provides the dispersion relation for slow, m = 1 
perturbations: 



S (1 — 2rf) uj — zu (uj + S) 



0. 



(26) 



whose solution is, 



-fc 



27j) ± ^\/S 2 (l - 2riY + 4tjj(uj 



S). 

(27) 

Equations (|26|l and (|27|) are invariant under (rj, u)) — > (1 — 
7/, —uj), because this operation is equivalent to interchanging 
the meaning of "prograde" and "retrograde" . It is convenient 
to first consider two special cases: 

(i) No counter-rotation, 77 = 0: When 77 = 0, equation (|27[1 
admits the two roots, uj = vj + S, an d uj = —w . The f ormer 
root corresponds to the eqn. (14) of iTrema inc (2001), and 
implies that the disc is locally stable to all m = 1 pertur- 
bations. However, uj — —uj is a spurious solution, arising 
from multiplication by in the derivation of equation (|18[) . 
Henceforth we assume that r\ 7^ 0. 

(ii) Equal counter-rotation, 77 = 1/2: When there is equal 
mass (locally) in the prograde and retrograde components, 
the two roots of equation (|27l) are, uj — ±\/uj (uj + S). If 
uj happens to be positive, then uj is real, and the disc is 
locally stable. However, uj < for most continuous discs, 
hence uj can be either real, or purely imaginary; there is no 
local overstability. The criteria for (in)stability are discussed 
below, along with the case of general 77. 

The sign of the discriminant of equation (|27p . 



V 



S 2 (l - 2n) 2 + 4uj(uj + S) 



(28) 



determines whether ui is real, or complex. If uj > 0, then 
T> > 0, hence uj is real. However, uj < for most continuous 
discs, and it is straightforward to determine that T> < 0, if 
5 lies in the range of values, < 5- < 5 < S+, where 

2\uj\ 



(i-2 V y 



[l ± 2^(1-7,) 



(29) 



However, we noted earlier that the maximum value that S 
can take is 5 max , given by equation (|23p . Hence, at a spec- 
ified r, T> is positive for all k, if S- > iS ma x. Therefore, the 
disc is stable to all short wavelength m = 1 perturbations, 
if 



61 



bo 



7rGS d 
ef2|ti7] 



(1 - 2t,) 2 



2\uj\ 



(1-2t?) 2 



n) 



Wvi^-v) 



(30) 



The ^-dependent factor in equation (|30p does not vanish for 
any value of 77. In fact, as n — > 1/2, the term in [ ] approaches 
unity, giving the stability criterion for the case 77 = 1/2, as 
may be verified independently. Since ([^[/fi) = 0(e) <S 1, 
the 6 required for local stability, according to equation (|30p . 
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equals bo (which is the minimum softening required for local 
axisymmetric stability) multiplied by a large factor, of order 
1/s. To the extent softening mimics "heat" (more precisely, 
the epicyclic radius) in collisionless discs, this criterion sug- 
gests that a Keplerian disc would have to very hot indeed, to 
be able to avoid a local instability to m — 1 perturbations. 
Hence one is led to consider overstable perturbations. 

Overstability occurs when T> < 0. We write cj = 
Q p ± iT, where Q p is the pattern speed of the m = 1 per- 
turbation, and F > is the growth rate. From equation (1271) . 



(1 - 2??) , 



§VI* 2 (i 



2 V y 



4ri7 (zu 



(31) 



(32) 



5 GLOBAL UNSTABLE MODES FOR EQUAL 
COUNTER-ROTATION 

We have seen in the previous section that in the case of 
equal counter-rotation (i.e. rj = 1/2) local analysis predicts 
purely unstable slow m = 1 modes. We will now go be- 
yond the local analysis and construct global m = 1 WKB 
modes for this case. To do this it is necessary to consider a 
concrete example. It is useful to take the axisymmetric un- 
perturbed disc to be a Kuzmin disc because (i) the Kuzmin 
disc is centrally concentrated and is hence a plausible candi- 
date for being a quite generic case; (ii) the surface density, 
the self-gravitational potential and the precession rate all 
have explicit analytical form s; (iii) the slow modes of the 
Kuzmin disc were studied in iTremainel (1200 ll ) for the case 
of no counter-rotation. 

The surface density of the Kuzmin disc is given by, 



Mr) 



aM, 



27r(r 2 +a 2 ) 5 / 2 ' 



(33) 



where M d is the disc mass and a is the central concentration 
parameter. The precession rate due to the Kuzmin disc is 
given by 



rod = 



3GM d a 2 



2fi(r)(r 2 + a 2 Y>/ 2 ' 



(34) 



where the rotational frequency is given by the Keplerian flow 
due to the central mass M as fi(r) = y GM/r 3 . Let us first 
consider the stability of the disc under axisymmetric pertur- 
bations. In section §[4] we derived the minimum value of the 
softening parameter bo that ensures local stability. In the 
slow mode limit k{t) = Vl[r) to the zeroth approximation. 
Therefore we find 



b 



2nGT, d aM d 



eM (r 2 + a 2 ) 3 / 2 



(35) 



The largest value of bo = aMd/eM, therefore the smallest 
softening parameter b that ensures stability everywhere sat- 
isfies b/a > Md/eMo- Let us define a parameter R = b/a 
and a parameter j3 through M d /Mo = PR. In terms of these 
parameters the previous inequality gives f3 < e. Substituting 
77 = 1/2 in equation (|27JI. we get 



5 = 



(36) 



Now we substitute for vj and recall that we are looking for 
global unstable modes. Therefore lu 2 = -r 2 , and 



a\k\ exp (— |fe|6) 



4/M0V r 2 (z* 2 + i) 4 



3 \MdJ fi(a) 



where x+ — r/a. By defining fc* 
V^T^r / '/3n(a) we finally obtain 



3 1 +x+ 2 ' 

(37) 

a\k\ and r* = 



fl|fc*|exp(-|fe*|.R) 

*(x* ; r*,fl) 



r*, R) , where 

r^ 2 (x* 2 + l) 4 3R 
R xj + l+a;* 2 



(38) 



5.1 Numerical Results 

Global unstable modes are determined by numerically solv- 
ing equation (|38p for fc* for a given value of F* and R and 
applying a quantization condition to obtain growth rate. 
Note that the right hand side of equation (|38[) blows up at 
x-t, = and x+ — > 00; however it is bounded from below 
and has a minimum at 1* ~ 1, the exact value depend- 
ing on the precise value of T* and R. Since the left hand 
side of this equation has a maximum value equal to unity 
at hi, — 1, the equation does not admit any solutions if 
min[$(a;*;r*,ii)] > 1. However, for min [4>(a;*; F*, ii)] < 1, 
the equation admits two real roots, one each on either side 
of ki, = 1. We denote the roots with fc* < 1 as the long- 
wavelength branch and the one with fc* > 1 as the short- 
wavelength branch. 

As noted above, the right hand side of equation (|38|) 
is unbounded from above and blows up at small and large 
values of x t . It is clear that the real roots exist only for a 
finite range of the radial coordinate, a < x+ < b, where both 
a and b depend on the parameters F* and R in a compli- 
cated manner. We shall assume that at these points the wave 
is reflected, and we therefore impose the Bohr-Sommerfeld 
quanitization condition 



f 



k*dx t 



n + 



1 



(39) 



We find that the quantization condition cannot be satis- 
fied for R > 1. However, R = b/a > M d /eM 1; therefore 
a small value of R is allowed by the inequality. We seek the 
global modes for two cases; R — 0.1 and R — 0.01. In Fig- 
ure |T]) and J3} we display the phase plots for the short and 
long wavelength branches for R = 0.01 and T* = 0.0001. 
The general trend is that as T* increases the contours be- 
come narrower in the horizontal direction and expand in 
the vertical direction in such a manner that the area un- 
der the closed curve decreases as F* increases. This can be 
seen clearly in Figures @ and @, where we have plotted 
the spectrum for the short and long wavelength branches, 
where we see that larger quantum numbers correspond to 
smaller values of the growth rate. The two figures show that 
the effect of changing R is to translate the spectrum hor- 
izontally. This is only an approximate behaviour and does 
not imply anything special. More signicantly the spectrum 
is seen to be very close to a power law. This behaviour is 
quite robust and persists for other values of R. 
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Radial coordinate x t Radial coordinate x„ 

Figure 1. Phase plot of a possible mode in the short— wavelength Figure 3. Phase plot of a possible mode in the long-wavelength 

branch for R = 0.01 and I\ = 0.0001. branch for R = 0.01 and I\ = 0.0001. 
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01 




Quantum Number n 

Figure 2. Dimcnsionlcss growth rate versus quantum number of 
the mode for the short— wavelength branch. 



6 CONCLUSIONS 

The principal aim of this work is to present the simplest 
model that permits a largely analytical exploration of the 
m = 1 counter-rotating instability in a "hot" nearly Keple- 
rian disc of collisionless self-gravitating matter. To this end 
we have considered a two-component softened gravity disc, 
and performed a linearised WKB analysis of both local and 
global modes. We derive an analytical expression for local 
WKB waves for arbitrary m, which turns out to be quartic 
in the frequency u). Specialising to m = 1, we show that 
uj is smaller than the (Keplerian) orbital frequency by the 
small quantitity e = Md/M (the ratio of the disc mass to 
the mass of the central object); in other words, the m = 1 
modes are slow modes. The dispersion relation now reduces 
to a quadratic equation in lu. Hence the criteria for stability, 
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Figure 4. Dimcnsionlcss growth rate versus quantum number of 
the mode for the long— wavelength branch. 



instability and overstability can be readily derived in sim- 
ple analytical forms. For a one-component disc (which does 
not have any counter-rotation), t he m = 1 modes are sta- 
ble, consistent with the results of iTremaind (|200lh . Equal 
mass in the two counter-rotating components corresponds 
to the case of not net rotation. In this case we find that 
the local mode s are purely unstable (i.e. not over s table) , 
consistent with Arakil Jl987l ); Palmer fc Papaloizoul dl990T): 
ISellwood fc MerrittJ <|l994D ; iLovelace et al. I (ll997f ); Irbumal 
(|2002l ). However the general case of arbitrary mass ratio in 
the two counter-rotating components corresponds to over- 
stability, and we show analytically that the discs must be 
unrealistically hot to avoid an overstability. We finally con- 
tracted global WKB modes, numerically, for the case of a 
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Kuzmin disc for the case of no net rotation, by using Bohr- 
Sommerfeld quantisation. 
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